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Abstract 

The action of a Backlund-Darboux transformation on a spectral 
problem associated with a known integrable system can define a new 
discrete spectral problem. In this paper, we interpret a slightly gener¬ 
alized version of the binary Backlund-Darboux (or Zakharov-Shabat 
dressing) transformation for the nonlinear Schrodinger (NTS) hier¬ 
archy as a discrete spectral problem, wherein the two intermediate 
potentials appearing in the Darboux matrix are considered as a pair 
of new dependent variables. Then, we associate the discrete spectral 
problem with a suitable isospectral time-evolution equation, which 
forms the Lax-pair representation for a space-discrete NTS system. 
This formulation is valid for the most general case where the two de¬ 
pendent variables take values in (rectangular) matrices. In contrast to 
the matrix generalization of the Ablowitz-Ladik lattice, our discretiza¬ 
tion has a rational nonlinearity and admits a Hermitian conjugation 
reduction between the two dependent variables. Thus, a new proper 
space-discretization of the vector/matrix NLS equation is obtained; by 
changing the time part of the Lax pair, we also obtain an integrable 
space-discretization of the vector/matrix modified KdV (mKdV) equa¬ 
tion. Because Backlund-Darboux transformations are permutable, we 
can increase the number of discrete independent variables in a multi- 
dimensionally consistent way. By solving the consistency condition on 
the two-dimensional lattice, we obtain a new Yang-Baxter map of the 
NLS type, which can be considered as a fully discrete analog of the 
principal chiral model for projection matrices. 


1 


Contents 


1 Introduction 


a 


2 


A proper discretization of the matrix NLS hierarchy 

2.1 The continuous matrix NLS hierarchy. 

2.2 Backlund-Darboux transformation as a discrete spectral prob¬ 
lem . 

2.3 Isospectral time evolution . 


I 

[9 

[10 


3 Yang—Baxter map 


16 


4 Concluding remarks 


24 


A 


Integrable discretization of the vector sine-Gordon equation 


25 


References 


27 


2 











1 Introduction 


The problem of integrable discretization [T] has a close relationship with 
the theory of auto-Backlund transformations; this fact has been sporadically 
noticed for partial differential or differential-difference equations since the 
mid-1970s. Calogero and Degasperis w and Chiu and Ladik |3] showed 
for some specihc examples that a class of auto-Backlund transformations can 
be identified as discrete-time flows that belong to the same integrable hier¬ 
archy as the original continuous-time flows. Hirota [5] and Orfanidis EE] 
showed for the sine-Gordon equation in light-cone coordinates that a one- 
parameter auto-Backlund transformation and the associated nonlinear su¬ 
perposition principle (Bianchi’s permutability theorem) lead directly to its 
integrable discretizations. These results have culminated in the beautiful 
notion of Miwa variables (or Miwa shifts) in Sato theory [8l|9]. In rough 
terms, a one-parameter auto-Backlund transformation, called an elementary 
Backlund transformation, can be reinterpreted as a discrete-time flow, which 
can generate an inhnite number of continuous-time flows through the Taylor 
series expansion in the Backlund (or time-step) parameter; the commuta¬ 
tivity of two elementary Backlund transformations with (generally) different 
values of the Backlund parameter provides a fully discrete integrable system 
that admits a zero-curvature (or Lax-pair cni) representation and contains 
information on the continuous-time flows. Such a unified point of view is 
so fascinating that one is tempted to believe that proper discretizations of a 
given integrable system can always be obtained from its elementary Backlund 
transformation or the associated nonlinear superposition formula. However, 
in fact, it is not that simple; we remark the following points. 

• An integrable system appears as a member of an infinite hierarchy of 
commuting flows, which is actually a bi-infinite hierarchy comprising 
the positive and negative flows. For example, the sine-Gordon equa¬ 
tion is the first negative flow of the integrable hierarchy, while its first 
nontrivial positive flow is (the potential form of) the modified KdV 
(mKdV) equation [II11I2]. It is not evident which particular flows in 
a bi-in£nite hierarchy can be discretized by considering an elementary 
Backlund transformation and how to take the continuous limit. 

• The idea to interpret an elementary Backlund transformation and the 
associated nonlinear superposition formula as defining discrete inte¬ 
grable systems is useful for integrable hierarchies with one scalar un¬ 
known. However, it is not the case for two (or more) component systems 
that admit a complex conjugation reduction between a pair of depen¬ 
dent variables, such as the nonlinear Schrodinger (NLS) system [miT^ . 
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as well as their matrix generalizations that admit a Hermitian conjuga¬ 
tion reduction between a pair of matrix dependent variables, such as the 
matrix NLS system [13]. Indeed, an elementary Backlund transforma¬ 
tion generally does not maintain the complex/Hermitian conjugation 
reduction; thus, for a given integrable equation obtained as the result of 
the reduction, we cannot derive its proper discretization directly from 
an elementary Backlund transformation for the original nonreduced 
system (cf. [TH - [T7] j. Actually, we can consider a suitable composition 
of two (or more) elementary Backlund transformations [laiiHiiini so 
that the composite {e.g., binary) Backlund transformation can main¬ 
tain the complex/Hermitian conjugation reduction; however, such a 
composite transformation involves either a square-root function with 
an indehnite sign [201121] or a nonlocal operation such as an indehnite 
integral [3], so it does not lead to discrete integrable systems that can 
be written in local form and dehne a unique time evolution (cf. (7.18) 
in [22] and (4.17) in [23]). 

Thus, for each integrable system that admits the complex/Hermitian conju¬ 
gation reduction, we need to construct its proper discretization on a case-by- 
case consideration and accumulate more knowledge on this subject. 

In this paper, we construct a new proper space-discretization of the (gen¬ 
erally rectangular) matrix generalization of the NLS system, namely, the 
matrix NLS system [13], which is an integrable system admitting the Hermi¬ 
tian conjugation reduction between the two matrix dependent variables and 
contains the vector NLS equation known as the Manakov model [2l] as a 
special case. For the scalar NLS system, the proper and most natural space- 
discretization was found by Ablowitz and Ladik [251|2S] in the mid-1970s; 
however, its straightforward matrix generalization [27] (also see [28l|29] and 
references therein) does not admit a Hermitian conjugation reduction be¬ 
tween the two matrix variables in local form, so it is not a proper space- 
discretization of the matrix NLS system. To obtain a proper space-discrete 
matrix NLS system, we follow the new approach introduced in our previous 
paper That is, we hrst reinterpret (a slightly generalized version of) 

a binary Backlund-Darboux transformation for the continuous matrix NLS 
hierarchy as a discrete spectral problem [16], where the two intermediate po¬ 
tentials appearing in the binary Backlund-Darboux transformation are con¬ 
sidered as new dependent variables. Then, in view of the peculiar structure of 
the discrete spectral problem, we associate it with a suitable isospectral time- 
evolution equation to compose a Lax pair and derive an evolutionary lattice 
system from the compatibility condition called the zero-curvature equation. 
This lattice system involves some free parameters and, with a suitable choice 
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of the parameters, it admits the Hermitian conjugation reduction between 
the two dependent variables and provides a proper space-discretization of the 
matrix NLS system; thus, it can generate proper space-discrete analogs of 
various multicomponent NLS equations [T^I24[l31fl36] obtained as reductions 
of the matrix NLS system. 

The space-discrete NLS system derived in this paper has a rational non¬ 
linearity, so it is intrinsically more complicated than the Ablowitz-Ladik 
discretization of the NLS system |25l|26] that has a polynomial nonlinearity. 
However, the rational nonlinearity can exhibit a saturation effect, so in some 
sense it is more “physical” than the polynomial nonlinearity. 

The space-discrete matrix NLS system appears as a member of an inhnite 
hierarchy of commuting flows. Other flows of this integrable hierarchy can 
be obtained by changing the temporal part of the Lax pair. In particular, we 
can derive a space-discrete analog of the matrix mKdV system, which admits 
the Hermitian conjugation (or matrix transposition) reduction between the 
two dependent variables. Thus, integrable space-discretizations of various 
multicomponent mKdV equations |191I231I31[[3B] - I3^ can be obtained through 
reductions. 

We derive the space-discrete matrix NLS hierarchy from a (slightly gener¬ 
alized) binary Backlund-Darboux transformation for the continuous matrix 
NLS hierarchy by utilizing it as the underlying discrete spectral problem. 
Actually, we can consider an arbitrary number of Backlund-Darboux trans¬ 
formations with (generally) different values of the Backlund parameters and 
assign a new discrete independent variable to each Backlund-Darboux trans¬ 
formation; that is, a single application of each transformation is identified 
with a unit shift of the corresponding discrete independent variable. The 
consistency of such a multidimensional lattice is guaranteed by Bianchi’s 
permutability theorem for two Backlund-Darboux transformations on each 
quadrilateral. In our approach, the dependent variables are the intermediate 
potentials that appear in each Backlund-Darboux transformation, so they 
are assigned to the edges of the lattice [401144] . Then, the permutability 
condition for Backlund-Darboux transformations on a quadrilateral can be 
solved explicitly, providing a new Yang-Baxter map (to use Veselov’s termi¬ 
nology [inillS]); this is apparently related to the factorizability property of 
an Y-soliton collision into pairwise collisions in the continuous matrix NLS 
hierarchy (cf. [46]). The permutability condition can be regarded as a matrix 
re-factorization problem, which directly provides the Lax representation for 
the Yang-Baxter map [401144114411471 . We can identify this Yang-Baxter map 
as a fully discrete analog of the principal chiral model [4S] restricted to the 
space of projection matrices. 
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This paper is organized as follows. In section 2, we introdnce (a slightly 
generalized version of) the binary Backlnnd-Darbonx transformation for the 
continnous matrix NTS hierarchy, consider it as a discrete spectral problem 
and associate it with a snitable isospectral time-evolntion eqnation. Then, 
the compatibility condition for this Lax pair provides a lattice system involv¬ 
ing some arbitrary parameters; with a snitable choice of the parameters, it 
provides a proper space-discrete analog of the matrix NLS system. By chang¬ 
ing the choice of the parameters appropriately, we also obtain a proper space- 
discretization of the matrix mKdV system. In section 3, we solve Bianchi’s 
permutability condition for the Backlnnd-Darbonx transformations in such 
a way that a new Yang-Baxter map of the NLS type together with its Lax 
representation can be obtained; this considerably generalizes an analogous 
result of Goncharenko and Veselov [131IT7] . Section 4 is devoted to concluding 


remarks. 


In appendix we present some result obtained by the author in the 
late 1990s but left unpublished. We consider a Lax-pair representation for 
a vector analog of the sine-Gordon equation [IHIED] and describe how to 
discretize one of the two independent variables preserving the integrability; 
this is somewhat related to sections 2 and 3 of this paper. 

2 A proper discretization of the matrix NLS 
hierarchy 

In this section, we utilize a (slightly generalized) binary Backlnnd-Darbonx 
transformation for the continuous matrix NLS hierarchy as a discrete spectral 
problem and derive proper space-discrete analogs of the matrix NLS system 
and the matrix mKdV system. 

2.1 The continuous matrix NLS hierarchy 

In 1974, Zakharov and Shabat [13] proposed the matrix generalization of the 
NLS system [TT|P:2]: 


iQt2 + Qxx ~ ^QRQ — O, 

iRt2 ~ Rxx + "^RQR = O, 


( 2 . 1 ) 


where the dependent variables Q and R are li x I 2 and I 2 x li (generally 
rectangular) matrices and the subscripts denote the partial differentiation. 
The symbol O on the right-hand side of the equations is used to stress that 
the dependent variables can take values in matrices. The matrix NLS system 
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m, as well as related equations, has been studied intensively in recent years 
(see, e.g., [22] and references cited therein). 

The matrix NTS system fl2.ip is obtained as the compatibility condition 
for the overdetermined linear system of partial differential equations [5n[52] : 


Ti ■ 
T 2 

X 

■-iC4 Q 

R iC4 . 

' ■ 

^2 



Ti ■ 
^2 

t2 

— 2i((^^/q — iQR 2C,Q -\- iQx 

2(R — iRx 2i('^Ii^ -b iRQ 

1-1 

to 

1 _ 1 


Here, C is a constant spectral parameter, and J/j and are the h x li and 
I 2 X I 2 unit matrices, respectively. In the following, we suppress the indices 
of the unit matrices, noting that the dependent variables can take values 
in (generally) rectangular matrices. Equations 02.21) and 02.31) comprise the 
Lax-pair representation [101 for the matrix NLS system 02.11) . Any constant 
scalar matrix can be added to each Lax matrix, which does not affect the 
compatibility condition. 

The matrix NLS system 02. ip is a positive flow in a bi-infinite hierarchy 
of mutually commuting flows, which are all associated with the same spectral 
problem 02.21) . The next higher flow in the matrix NLS hierarchy is a matrix 
generalization [5T[[52] of the complex mKdV equation [TTlll21l53j . which reads 


Qti + Qxxx ~ SQxRQ — SQRQx — O, 
Rts + Rxxx ~ SRxQR — 3RQRx = O. 


The first negative flow in the matrix NLS hierarchy is a matrix analog of the 
complex sine-Gordon equation, which we rewrite in the formQ 

' Qt-i + 2n(/ + vu)~^ = O, 

Rt , — 2v(I + uv)~^ = O, 

) ^ > ( 2 . 5 ) 

Ux + 2iku = Q — uRu, 

Vx — 2ikv = —R + vQv. 

As long as we consider the first negative flow 02.5p for a single fixed value of k, 
the free parameter k is nonessential and can be set equal to zero by applying 

^ The first negative flow of the scalar NLS hierarchy mm admits a number of differ¬ 
ent expressions, which are related through simple transformations of dependent variables. 
Such different (but essentially equivalent) expressions are often called by different names 
in correspondence to different physical phenomena (or mathematical objects), such as self- 
induced transparency (or MaxwelLBloch) equations [T21[S3][^ (also see §4.4.b of [55]). 
stimulated Raman scattering, complex sine-Gordon equation, Pohlmeyer-Lund-Regge sys¬ 
tem, reduced nonlinear cr-model, etc. There are too many references to mention here. 
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a point transformation involving thns, it reduces to the simpler form: 


" <5t_i + 2n(/ + vu) ^ = O, 

Rt , — 2v(I + uv)~^ = O, 

^ ^ ( 2 . 6 ) 

Ux = Q — uRu, 

Vx = —R + vQv. 

However, it is often more convenient to leave k as a free parameter. Note that 
we can consider an arbitrary linear superposition of the hrst negative flow 
02.51) for different values of k [53] . In the case of scalar dependent variables, 
we can rewrite 02.51) (or 02.61) 1 as a closed non-evolutionary system for (Q, R) 
or (u,n). 

The time part of the Lax-pair representation 02.31) for 02.11) is now replaced 
with 


■ ■ 


^2 

tz 


—4iC^/ — 2iC,QR + QxR — QRx 
- 2iC,Rx - Rxx + 2RQR 


+ ‘2‘KQx ~ Qxx + 2QRQ 
-|- 2i(^RQ -|- RxQ — RQx 


for 02.4p and 


\ 1 

i 

r I 1 


\ -I u) 

-1 

r 1 

^2 

1 

1 

0 

+ 

V I 

) 

^2 


(2.7) 


for 02.5p . respectively. Here, the underscored part. 


P : = 


I 

O 

uv {I + uv)~^ u{I + vu)~^ 
v{I + uv)~^ {I + vu)~^ 


+ 


-I u 
V I 


u 

I 


(/ -f vu) ^ [ n / ] 


' -I u' 


■ 0 


' -I u' 

V I 


I 


V I 


( 2 . 8 ) 


is a projection matrix, i.e., it satishes the condition = P. Note, inciden¬ 
tally, that the free parameter k in 02.5p and 02.7p can be generalized to an 
arbitrary function of t_i. 

One of the most important properties of the matrix NLS hierarchy, which 
is crucial for its physical applications, is that it admits the Hermitian conju¬ 
gation reduction between the pair of dependent variables Q and R [33] (also 














































see pTl[52|IFrl - [5^ ): this reduction results in a self-focusing, self-defocusing or 
mixed focusing-defocusing nonlinearity, which is briefly summarized in §2.1 of 
our previous paper [60] . In particular, by imposing the reduction R = —Q^ 
on the matrix NLS system fl2.ip where the dagger denotes the Hermitian 
conjugation, we obtain the self-focusing matrix NLS equation |13j : 

^Qt 2 + Qxx + ‘^QQ^Q = O. (2-9) 

In the same manner, the reduction R = —Q^ simplifies fl2.4p to the matrix 
complex mKdV equation: 

Qti + Qxxx + ^QxQ^Q + ^QQ^Qx = O. ( 2 . 10 ) 

Another important property of the matrix NLS hierarchy is that the odd- 
order flows admit the matrix transposition reduction between the pair of 
variables Q and R. In particular, the matrix mKdV equation [231EZ113B], 

Qtz + Qxxx + “iQxQ^Q + “iQQ^Qx = O, (2-11) 

is obtained by imposing the reduction R = —Q^ on fl2.4p . where the super¬ 
script ^ denotes the transpose of a matrix. In the case of a vector dependent 
variable, fl2.1ip gives the vector mKdV equation [T91I3T] : 

+ Qxxx + 3(g, qQq + 3(q, q)q^ = 0. (2.12) 

If qr is a two-component real-valued vector, i.e., q = {qi,q 2 ), qi,q 2 ^ IR, the 
vector mKdV equation 02.121) can be rewritten as a single complex-valued 
mKdV equation for q := gi -|- ig 2 ; this is often referred to as the Sasa-Satsuma 
equation [6T] . 

Any proper discretization of the matrix NLS hierarchy should retain the 
feasibility of such reductions. In short, a discrete analog of the matrix NLS 
hierarchy is physically meaningful only if it is integrable and admits the Her¬ 
mitian conjugation (or matrix transposition) reduction similar to the contin¬ 
uous case. That is, it should be able to provide integrable discretizations for 
the reduced equations such as fl2.9l) - fl2.12l) . 


2.2 Backlund—Darboux transformation as a discrete 
spectral problem 


To obtain a proper discrete analog of the matrix NLS hierarchy, we start 
with a slightly generalized version of the binary Backlund-Darboux trans¬ 
formation [T6l[62l|63] for the continuous matrix NLS hierarchy as given by 


1 

1_ 


1- 

to 

"U 


{(5 + a)I 


<1 + 1^)1 


(ay -F ( 35 ) 


7/ 

u 

"1 

■ ■ 

V 

51 

/ 

^2 


(2.13) 
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Here, a, /?, 7 and 6 are arbitrary parameters, except that they are required 
to satisfy the condition a'y + f36 ^ 0 ; the special case 6 = — 7 ( 7 ^ 0 ) corre¬ 
sponds to the conventional binary Backlund-Darboux transformation [63] or 
the Zakharov-Shabat dressing method caEg, up to an inessential overall 
factor. Unlike the usual formulation of the binary Backlund-Darboux trans¬ 
formation, we do not express the intermediate potentials u and v in terms 
of linear eigenfunctions of the original Lax-pair representation [TUfTB] . Al¬ 
ternatively, we will consider u and n as a pair of new dependent variables 
dehned on a lattice, wherein the lattice index n G Z is understood intuitively 
as the number of iterations of the Backlund-Darboux transformation. To 
make it more explicit, we rewrite fl2.13p as a discrete spectral problem 


^l,n+l 

^ 2 ,n+l 


= Lr. 






2,n 


(2.14) 


where the Lax matrix is given by 


Lr^ 


{(6 + a)I 


(C<^ -|- Oi)I 


(-C7 + /5)/ 


- (ay -b (36) 


ll Un 

Vn 61 


-1 




-Cl + (3)I _ 


- (ay -b (36) 

6{'j6I - UnVn) ^ -{'y6I-UnVn) 

-(^61 - VnUn)~^Vn j(-/6I - 

61 Uyi 

_ -Vn 7 / 

1 -1 

--1 

+ 

T 

1 

v- 

1 _1 

y/ Un 
_ Vn 61 _ 

7 / Un 
Vn 61 

-1 

1 -1 

+ 

T 

1 

1 _1 

61 —Ur 

-Vn 'll 


-1 


(2.15) 


Actually, the Backlund parameters a, (3, 7 and 6 satisfying the condition 
ay -h (36 0 can be arbitrary functions of the discrete independent variable 

n, but for brevity we usually consider them as constants. Each of the four 
expressions in fl2.15p has its own advantages. 


2.3 Isospectral time evolution 


To compose a Lax pair, we need to associate fl2.14p with a suitable isospectral 
time-evolution equation. 


'hi,„ 


Mn 


^ 2 ,n 


(2.16) 
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Here, the spectral parameter ( involved in and is isospectral, i.e., 
Ci = 0. The compatibility condition for fl2.14p and 02.161) is given by (a space- 
discrete version of) the zero-curvature equation [25112611641165] 


Ln,t + LnMn — M^+lLn — O, 


(2.17) 


where and Mn comprise the Lax pair. For the Lax matrix in 02.151) . 
we use the hrst expression to compute with the aid of the formula 
= —X~^XtX~^ valid for a square matrix X so that the equations 
of motion for Un and Vn can be obtained explicitly. 

Let us hrst consider the original continuous spectral problem 02.2p . which 
now reads 

_ ^2,n 

The Backlund-Darboux transformation 02.141) with 02.151) preserves the form 
of the spectral problem 02.181) invariant. In the general case where the 
Backlund parameters a, (3, 7 and 6 depend on n (a ^ etc.), the com¬ 
patibility condition (i.e., 02.17p with t replaced by x) provides a system of 
differential-difference equations: 

'Jn-Qn+l T ^nQn T 2i(Q;^7n T l3n^n)'^n{Xn^nI 'VinUn) 

T T 2i(Q:jj7n T l3n^n)'^niS)n^n3^ ^n^n) O, 

(cUnyn T /3n^ri)'^n,x T Cyn'JnQn+l f^n^jiQn ’Unifi^riRn+l /3nRn)'^n 

T /3n^n}^n,x T ^n^l^nQn+l ^nQri)^n O. 

(2.19) 

When 7 ^ 0 , the hrst two equations relate the new potentials Qn+i and 
Rn+i to the old potentials Qn and Rn in the spectral problem 02.18p through 
the intermediate potentials Un and Vn (cf. [E]); in this case, 02.191) can be 
reformulated as Riccati equations for and 



-KI 

Qn 

^l.n 

X 

Rn 

KI 

^2,n 


(2.18) 


^nUn,x ‘Roin'^n T ^uRn^ri) 

'yn'^rijX T ^nQn'^m 

^nUn,x 2i/3fi'Ufi -|- 7nQn-|-l UnRn+lUni 
^ ^n^n,x 2iQ(jj'yjj -|- 6^Rn-\-l V^Qn+l^n- 


( 2 . 20 ) 


Thus, the intermediate potentials and can be expressed in terms of the 
linear eigenfunctions of the spectral problem 02.181) [lUIIS], e.g. 




-1 I 
n+1 /•— otn 
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where and ^i,n+i are square matrices. In the case of scalar dependent 
variables, we can eliminate Qn and Rn from fl2.20p to obtain a closed system 
for Un and 

With a suitable choice of the parameters, 02.191) can be considered as 
an integrable time-discretization of the matrix sine-Gordon system 02.61) . 
Indeed, in the case of /3„ = —and 6n = —7n, 02.19P provides a (generally 
nonuniform) discrete-time analog of 02.6p : 


Qn+l Qn -|- VnUn') G, 

Rn+1 Rn “ 1 “ ( 7 ,^/ -|- Un'Vn') G, 

< 2 (2-21) 

‘^'yn'^n,x 7n(Q^+l QQ ^ni^Rn+1 R Rn)^n G, 

^ ‘^'yn^n,x 7n('^n+l -^n) d* '^ni,Qn+l d“ QQ'^n G. 

To obtain proper space-discretizations of positive flows of the matrix NLS 
hierarchy, we consider a temporal Lax matrix that matches well with the 
factorized form of the spatial Lax matrix in 02.151) . Note that the hrst 

expression in 02.151) implies the C-independence of Ln,t, so LnMn — Mn+iLn in 

02.17P is also ^-independent. A simple ansatz for satisfying this condition 
automatically is 


= 


7^ Un 

Vn SI 


1 F 


C^-l-o 


G. 


7 / 

Un—l 


Un—1 

61 


+ 


cl 


dl 


( 2 . 22 ) 


where and G„ are C-independent square matrices and c and d are (- 
independent parameters. The (^-dependence of the main part of in 02.22p 
implies that it can be considered as a linear superposition of the temporal 
Lax matrix in 02.71) at two particular values of k (up to the addition of a 
constant scalar matrix), i.e., it corresponds to the two “first” negative flows 
of the matrix NLS hierarchy. Substituting 02.15P and 02.221) into 02.17p . we 
obtain recurrence relations for and G„; to satisfy them identically, we set 

Fn = a (07 -1- (35) (7^/ Un-lVn) \ Gn = b (07 + (36) ( 5 ^/ -h Vn-lUn) ^ , 


where a and b are arbitrary constants. In fact, a and b, as well as c and d in 
02.22p . can depend arbitrarily on the time variable t, but we do not consider 
it in this paper. With the above choice of Fn and Gn, 02.171) for 02.15p and 
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fl 2 . 22 p provides an evolutionary system of differential-difference equations: 


' Un,t + a - UnVn) I + Un-lVn) ^ {'JUn + 

-I- b {jUn+l + SUn) {S‘^1 + VnUn+l) ^ {'jSI - VnU^) “ (c - d)Un = O, 

< 

Vn,t - b {'j6I - VnUn) {5Vn + 7^^n-l) 

^ - a ((5n„+l + 7 n„) ( 7 ^/ -1- UnVn+l) ^ {ibl - UnVn) + (c - d)Vn = O, 

(2.23) 

which is independent of a and (3. Note that the order of products in fl2.23p 
can be changed using identities such as 


(7(5/ - UnVn) ( 7 ^-^ + Un-lVn) (7Wn + dUn-l) 

(5^ 1 1 

= 6Un H- Un-1 - {jUn + SUn-l) Vn I + (7«n + SUn-l) 

7 7 

= {jUn + SUn-l) + VnUn-l) ^ {'jSI - VnUn) ■ 


In the case of scalar Un and fl2.23p is related to known integrable 
systems such as the lattice Heisenberg ferromagnet model [661 - I69] through 
simple changes of dependent variables (also see (2.19)-(2.20) in [TD] together 
with im ); thus, it is not really a new integrable system. However, fl2.23p 
in the general matrix case, its reductions and identification as a discrete 
analog of the matrix NLS (or mKdV) equation have not been reported in the 
literature. 

Considering its relation to the space-discrete Kaup-Newell system (see 
Propositions 2.1 and 2.2 in [72] with a variable change —)■ v~^), we can 
show that the lattice system fl2.23|) possesses an ultralocal (but noncanonical) 
Poisson bracket. For scalar Un and it can be written as Un,t = {un, H} 
and Vn,t = {vn-, //} with the Hamiltonian and the Poisson bracket given by 


n 


a log 


/ + Un-lVn \ 

V 'yd -UnVn ) 


+ b\og 


/ 6^ +Un+lVn \ 

V 7 ^ - UnVn J 


+ 


c — d 
'yd - UnV. 


and 

\Unii Un} }Vmj Vn} 0 , {Umi Vn} dnin ( 7 *^ VnVn) 5 

respectively. Here, d^n is the Kronecker delta, which has nothing to do with 
the free parameter d. 

By rescaling the variables and parameters as Un —)■ 'yun, Vn —)■ 7 n„, a —)■ 7 a, 
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V 


7 


(2.24) 


Then, by taking the limit 7 —?• 0, 02.241) rednces to a known integrable space- 
discretization of the NLS system [271173]: 



T O , 


^ - dVn = O, 


'^n,t r (^n—1 ^n^n^n) nn+l (d “F 


which is associated with an elementary auto-Backlund transformation for the 
continuous matrix NLS hierarchy [HHIS] as outlined in na. This limiting 
case no longer admits the complex (or Hermitian) conjugation reduction be¬ 
tween Un and Vn, so it does not provide a physically meaningful discretization 
of the NLS system. 

When the parameters satisfy the conditions h = 7 *, 6 = —a* and d = c*, 
the lattice system fl2.23p admits not only the complex conjugation reduction 
for square matrices Un and Vn but also the Hermitian conjugation reduc¬ 
tion for (generally) rectangular matrices Un and In particular, in the 
case of h = 7 = 1, 6 = —a* and d = c*, the Hermitian conjugation reduction 
Vn = simplihes 02.231) to 

® (d ^n+l^n) (^n+l T ^n) 



O. (2.25) 


By further setting a = —i and c = —2i, 02.251) provides a new integrable 
space-discretization of the matrix NLS equation 02.91) . 



(2.26) 
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This equation is defined on the three lattice sites n — 1, n, n + 1, so it is 
simpler and appears to be more interesting than the space-discrete matrix 
NTS equation proposed in our previous paper [28] (see (4.1) therein), which 
depends on five lattice sites after imposing the Hermitian conjugation re¬ 
duction. Considering further reductions of the matrix dependent variable 
Un, we can obtain integrable space-discretizations of various multicompo¬ 
nent NTS equations |31H3S]- In particular, when Un is a row (or column) 
vector, 02.261) provides a new integrable space-discretization of the vector 
NLS equation often referred to as the Manakov model [23|. Equation 02.26|) 
corresponds to the self-focusing case, but it is also possible to consider the 
self-defocusing case or a mixed focusing-defocusing case as in the continuous 


Note that a fully discrete matrix NLS equation was proposed by van der 
Linden, Nijhoff, Capel and Quispel in 1986 (see (4.17) in [23]). However, it 
involves a square-root function of a general non-diagonal matrix defined using 
the Taylor series expansion, so it is not an explicit and closed expression. In 
fact, even in the scalar case, their fully discrete NLS equation (see (7.18) 
in [2^) is an implicit equation, which does not define the time evolution 
uniquely. 

In the case of a = — 1 and c = 0, fl2.25|) provides an integrable space- 
discretization of the matrix complex mKdV equation 02.101) given by 

'^n,t T T (-f ('^n+l T ^n) 

“ 1 “ Unul^ (7 Un—lU^n) i.^n T l) O. 

When the parameters satisfy the conditions 6 = 'y, b = —a and c = d = 0, 
the lattice system 02.231) admits the matrix transposition reduction between 
Un and Vn- In particular, by setting 5 = 7 = 1, b = —a = 1, c = d = 0 and 
Vn = 02.231) reduces to an integrable space-discretization of the matrix 

mKdV equation 02.111) : 

Un,t T (7 -|- UnUn) (7 Un+lUn) (^Un+1 T Un) 

(7 -|- UnUn') (7 Un—lUn') {Un T l) O. 

In the vector case, this reads 


'^n,t T (1 -|- (yu,ni Un)) 


'^n+1 T 

1 '^n) 


Un Un—1 
1 '^n—l) 


0, (2.27) 


which provides an integrable space-discrete analog of the vector mKdV equa¬ 
tion 02.12p . In the single-component (i.e., scalar) case, this can be rewritten 
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in a more familiar form (cf. (3.5) in [75]) by setting Un ='■ tanra^. When 
Un is a real two-component vector, i.e., = {un\un'’)i Un\un'’ € K, we 

can introduce a complex-valued function '■= Un'^ -|- lUn^ to rewrite fl2.27p 
as a single evolution equation for Then, we obtain an integrable space- 
discretization of the complex mKdV equation often referred to as the Sasa- 
Satsuma equation [ 6 T| . 

When 6 = ± 7 , b = —a and c = d = 0, we can reduce the lattice system 
fl2.23p to a single closed equation for by setting Vn as a constant scalar 
matrix. In particular, by setting S = b = —a, c = d = 0, Vn = — 7 ^/ and 
07 = —1, we obtain a proper space-discretization of the matrix KdV equa¬ 
tion [10] (cf. fl2.4p with R = const.), which reads 

^n,i T (^n+l T (^n+l T '^n) (7 Un+l) (^n+l T '^n) 

T ^n—l) (7 M^—i) {Un T l) O. 

Alternatively, by setting 5 = — 7 , b = —a, c = d = 0, = 7^7 and 07 = 1, 

we obtain 


Un,t T (^n+l T Un—\ 2'Uji) (Un-\-\ Un) (7 -|- Un-\-l) (^n+l Un) 

(Un Un—l) (7 -|- Un—l) {Un l) 

which is an integrable space-discretization of the equation 

Ut Uxx ^Ux^d u) Ux O 

that can be linearized as 

[(/ + «)-'],+[(/ + «)-■]„ = 0 . 

3 Yang—Baxter map 

In this section, we utilize the permutability property of Backlund-Darboux 
transformations to construct a new Yang-Baxter map, which provides a fully 
discrete analog of the principal chiral model for projection matrices and ad¬ 
mits the Hermitian conjugation reduction. 

Let us consider two copies of the “same” Backlund-Darboux transforma¬ 
tion with different sets of Backlund parameters and intermediate potentials 
and assign them to two lattice directions m and n, respectively. A forward 
shift in the n-direction is defined as in fl2.14p with fl2.15p . i.e. 

^ m,n+l m,n (3.1) 


16 











with the Lax matrix Lm,n given by 


L 


m^n 


(C<^ + tt)/ 

(-C7 + /5)/ 


(ay + (35) 


yL Um,n 

Vm,n 51 


(3.2) 


Here, the condition ay + 7 ^ 0 is assnmed; for brevity, the n-dependence 

of the Backlund parameters a, (3, ■y and 5 is snppressed. We are only inter¬ 
ested in the general case of y5 7 ^ 0 and do not consider the limiting case of 
y (or 5) —)■ 0 where the binary Backlnnd-Darboux transformation can rednce 
to an elementary Backlund-Darbonx transformation. 

Similarly, we define a forward shift in the m-direction as 




m+l,n 


= K 


m.n^ m.n 


(3.3) 


with the Lax matrix Vm,n given by 


^m.n. — 


a)/ 

{-Cc + b)I 


(ac + hd) 


cl 

'^m,n 


Qm,n 

dl 


(3.4) 


Note that the parameters a, 6 , c and d used in this section have no direct 
relationship with those used in subsection 12.31 The condition ac -|- fed 7 ^ 0 is 
assumed and the m-dependence of a, b, c and d is suppressed. In addition, 
we assume cd ^ 0 and do not consider the limiting case of c (or d) —)■ 0 where 
the binary Backlnnd-Darboux transformation can reduce to an elementary 
Backlnnd-Darboux transformation. Note that the Lax matrices in fl3.2p and 
fl3.4p are both (/i -|- I 2 ) x (/i -|- I 2 ) block matrices, i.e., they are partitioned in 
the same manner. 

The compatibility condition for (13.1 p and fl3.3p is given by (a fully discrete 
version of) the zero-curvature equation [bl fTHl 1^ 176] 


^m+l,n ^m,n 




(3.5) 


The Lax matrices Lm,n and 14rx,n are of the same form (cf. (13.21) and (13.41) ). 
so this is a matrix re-factorization problem [1SIS3117]. This re-factorization 
problem can be solved explicitly; that is, for the Lax matrices originating 
from a Backlnnd-Darboux transformation as given above, the two matrices 
on the left-hand side of (13.51) can uniquely determine the two matrices on 
the right-hand side and vice versa. Thus, this provides a Yang-Baxter map 
admitting the Lax (or zero-curvature) representation [40114311441147] . 

With the aid of a gauge transformation 




m,n — 


5^d^I 

(-y)"(-c)'"J 


T' 

m,n’ 


(3.6) 
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we can reduce the general case of 7(5 7 ^ 0 and cd ^ 0 to the simpler case 
of (5 = —7 = 1 and d = —c = 1 after a minor redehnition of the parameters 
and the dependent variables. The Lax matrices in this case take the same 
form as the standard binary Backlund-Darboux transformation [63] or the 
Zakharov-Shabat dressing method [I6l[62| up to an inessential overall factor, 
so they admit a compact expression in terms of a projection matrix. 

Indeed, up to an inessential m-independent overall factor (cf. 03.51) 1. we 
can express the Lax matrix Lm,n in the case of (5 = —7 = 1 as 


Lm,n 


(C + 


(C + /9)/ 


{oi — /5) 


-/ Ur 


^Tn,n 

^m..n. I 


-1 




(3.7) 


Here, the projection matrix Pm,n is dehned as 


P ■= 

m..n. • 


O 


+ 


-I 

'^m.n 


U 


n -1 


m,n 

I 




^m,n 

I 

{I + v 

m,n^m,n) 

Vm^n 

l] 

-I 

^m,n 

' 0 

-I 

^m,n 

'^m,n 

I 

I 

'^m,n 

I 


(3.8) 


which indeed satishes {Pm,nY = Pm,n- Note that the correspondence between 
Pm,n and (um,n,Vm,n) in 03.8^ is oue-to-oue. Owing to the relation 



a —13 

C+^ 




13 — a 

C + 





the inverse of the Lax matrix Lm,n has the same form as Lm,n with a 77 /5, up 
to an inessential overall factor. Note, incidentally, that [±(/ — 2Pm,n)]^ = P 
Thus, using the “spin matrix” Sm,n ■= —I + 2Pm,n) the Lax matrix Lm,n can 
be rewritten in the form: 

Pm,n X / T —— ■ ■ '^Smm {,Sm.,n) P 

2Q + a + /3 

In the simplest 2x2 case, this is essentially the spatial Lax matrix for the 
lattice Heisenberg ferromagnet model [6611^7] (see [69] and §3.4 of [77] for the 
matrix generalization). 
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Similarly, we can express the Lax matrix Vm,n in the case of c? = —c = 1 


as 


^m,n 


(C + 


(C + h)I 


+ (a - 6) 


I Qm,r 

^ m.r). I 


-1 


0. — b ^ 
oc / + — 

C + b 


(3.9) 


Here, the projection matrix ^rn,n is defined as 


m,n • 


o + 

I Qm,i 
m,n I 


I Qm,n 

r* m..n. I 


-1 



■ 0 


I Qrn^n 


I 


f'm,n I 


-1 


(3.10) 


which satisfies {^rn,nY = ^m,n- 

Snbstitnting 03.71) and 03.9p into the fnlly discrete zero-curvatnre equation 
03.51) . we have 


T ^^ p 


or equivalently, 
^ b — a ^ 

I + 


/ + 


a — b 


, ^ a — b ^ 

= ( -f + ^ ^ ^m,n+l 


C + a 


/ + 


(3 — a 
C + a 


P, 


m+l,n 


= / + 


(3 — a 
C + a 


i + ^^p 

C + /3 " 

(3.11) 


^ b — a ^ 

^ ~Z ' I 7 

C + a ' 

(3.12) 


which generalizes the re-factorization problem studied by Goncharenko and 
Veselov [131H7] . This is an identity in the spectral parameter (, so at 0(1/C) 
it provides the conservation law: 


(ex /3) (Pm+l,n Pm,n) (^ b'j ^m,ri) ■ 

Using this conservation law, we can eliminate one of Pm+i,n) Pm,n: ^m,n+i 
and ^rn,n in the matrix re-factorization problem 03.lip (or 03.12p L Then, 
it can be solved explicitly so that the map {Pm,n, ^m,n+i) t {Pm+i,n, 3P^m,n) 
and the map {Pm,n, 33^m,n) t {Pm+i,m 33^m,n+i) as well as their inverses can be 
expressed in closed form; such a map represents a similarity transformation 
that transforms a projection matrix into a projection matrix of the same 
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rank. For instance, the map {Pm,n, ^m,n+i) ^ {Pm+i,n, ^m,n) reads 


Pm+l,n — I I 


a — (3 a — h 


a — b 


P 


a — b 
a — b 




m,n+l 1 ^m,n 


P I 

m..n. i 


a — (3 a — b 


a — b 


P 




^ = I / - - __ 

^ m.n 1 7 m,n 7 

' a — 0 a — 0 


t-x r) 


m,n+l I m,n+l 


a — b 
(3.13a) 

a — b 


m,n+l 


T _ ^^ p _ Z_1 

a — 0 a — 0 

(3.13b) 


which is a new Yang-Baxter map acting on the space of projection matrices 
and admits the Lax representation fl3.12p . If 03.13p is viewed as a system 
dehned on the two-dimensional lattice, the parameters a and [3 (resp. a and 
b) can be arbitrary fnnctions of the discrete independent variable n (resp. m). 
Thus, 03.13|) provides a new fully discrete analog of (a nonautonomous exten¬ 
sion of) the principal chiral model |1H] restricted to the space of projection 
matrices (P^ = P, = <^): 


dP 

dr] 




fiv) 


«(0 - a{v) 
9(0 

«(0 - a{v) 


[P, , 

[P, . 


(3.14a) 

(3.14b) 


Here, [P, := P3^ — 3^P is the commutator. When jij]) and are 

purely imaginary and (y{^) and 0 ( 77 ) are real, the principal chiral model fl3.14p 
admits the Hermitian conjugation reduction Pl = P, = 3^. 

The Lax-pair representation for 03.141) is 




C - a(0 


PT, 


^T] 


fiv) 

C - a{v) 




(3.15) 


where C is a constant spectral parameter. Note that we can normalize g{^) 
and f{g) to nonzero constants using a point transformation of the form 
X(,^), g ^Y{rj). This Lax-pair representation can be identified with 
02 .7p for two different values of A;, so we can understand the principal chiral 
model 03.141) as the compatibility condition for two “first” negative flows of 
the matrix NLS hierarchy. 

The nonautonomous principal chiral model 03.14p implies the relation 


I dP I 

fiv) dg g{0 


(3.16) 
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Thus, substituting the expression (cf. the Lax-pair representation fl3.15|) i 


P = 


a{7]) 


9(0 ’ fiv) 

into fl3.16p . we obtain the simpler form: 


4 , 4 -' 


a(() (454-‘),^ - o(r,) (4,4-‘)^ = O. 

It remains an open question whether this model is related to a more familiar 
nonautonomous chiral model that appears in general relativity (see, e.g., [TS] 
and references therein). 

Let us derive a more explicit component form of the Yang-Baxter map 
fl3.13p . In view of the representations of the projection matrices in fl3.8p and 
fl3.10p . we can obtain from fl3.13p the following relations: 



-I 








Y 

m,n 




J 

a — {3 


(/ 

H“ ^7n,n^m,n) 

I — 

a — b 

Vm 

,n “1“ 


)"^ (/ 

+ 


Qm 

n+l^m,n+l “1“ Q^m,n+l^m,n+l) 

Qm,n-\-l 

+ 

^m,n+l9m,n+l ) 


“1“ Qm 

,n+l^m,n+l) 

— i 

{I + r 

m, 

a+l9m,n+l') 

I 9m,n 

r 7 





Tm 

n I \ 

_ 

Wm,n _ 





J 

a — (3 


(I 

^7n,n^m,n) 

1 — 

a — b 

Vm 

,n “t" 

'^m,n 

Y (/ 

+ 

\-l 

Qm 

n+l^m,n+l “1“ Qm,n-\-l'^m,n 

+ 1 ) 


+ 

^ m,n-\-lQm,n-\-l) 


Qm 

,n+l^ m,n+l) 

— 1 

{I + r 

m, 

ri+l9m,n+l') 


a — b 


■X 


~'m,n 


a — b 

-rx 

a — b 
-I 




where Ym,n, Zm,n and Wm,n are some square matrices. By eliminating 

these unknown quantities, we can express the Yang-Baxter map in compo- 


'm,n 

I 


9m,n+l 

I 
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nent form as 


^m+ljn 




X 


I - 


a — b 
f3-b 

a — b 
(3-b 


Q.m^n+1 {I + rm,n+iqm ,n+l 




-1 


[(/^ ^m,n (-^ “l~ '^m,n+l^m,n) Q^m,n+1 (-^ ^m,nH-lQ^m,nH-l) ] 

^ [(/3 (-f H“ ^m,n+l^m,n) (-^ H” ^m,n+l?m,n+l) ] 

[(/3 (-f “1“ '^m,n^m,nH-l) &) (-f H“ Q^m,n+l^m,nH-l) ] 

X [(/3-6)(/ + (a-6) (J + 

Q'm,n+l^m,n+l) *?m,n+l] ) 

^m+l,n ®) '^m,n (-^ “1“ Qm,n+l^m,n') (® ^m,n+l (-^ “1“ ?m,n+l^m,n+l) ] 

^ [(Q^ ®) (-^ “1“ *?m,n+l^m,n) (® (-^ “f* Q'm,n+l^m,n+l) ] j 

qm,n ®) ?m,n+l (-^ “1“ ^m,n?m,n+l) /3) '^m,n (-^ “1“ '^m,n^m,n) ] 

X 1^(0 o) (/ + 'l’m,n<?m,n+l) (c^ / 3 ) (/ -|- Vm,nUm,n) ] 

[(® ®) ('^ Qm,n+l^m,n) (® /3) (-^ “1“ ] 

X 1^(0 o) (/ -|- qm,n+l'^m,n) qm,n+l (o^ /3) (/ ~l“ Um,n^m,n} ^m,n] ) 

f'm,n ~ [(/3 b) Tfn^n+1 (-^ ~l“ Um,nXm,n+l^ "i" (ci /3) Vm,n (-^ "1“ Um,nVm,n) ] 

^ [(/3 b') (/ -|- Urn,n^m,n+l') (o^ /3) (/ “1“ ] 


Note that each of the expressions for Um+i,n, Vm+i,n, qm,n and can be 
rewritten in many different (but nontrivially equivalent) forms, as is illus¬ 
trated for Um+i,n- In fact, these expressions can be derived more directly 
from the definition of the standard binary Backlund-Darboux transforma¬ 
tion [T6ll63] as well as its inverse, because the intermediate potentials in each 
projection matrix {e.g.^ Um,n and Vm,n in 03.81) ) can be expressed explicitly in 
terms of the linear eigenfunctions. In short, a companion map (see mmi!) 
represents the composition law for two binary Backlund-Darboux transfor¬ 
mations. 

When (3 = a* and b = a*, we can impose either the complex conjugation 
reduction Vm,n = rm,n+l = (rq*m,n+l, Wl,n = ^m.n = (^q*m,n 

for square matrices or the Hermitian conjugation reduction (* —)■ f) for (gen¬ 
erally) rectangular matrices, where a is a nonzero real constant; for simplic¬ 
ity, we do not discuss a more general reduction involving constant Hermitian 
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matrices corresponding to the mixed focusing-defocnsing case (cf. [niH33l[S3" 
EH]). In particnlar, by setting /3 = a*, b = a*, Vm,n = ’’m.n+i = qL,n+i^ 

Vm+i,n = '^rn,n = q}n,ny obtain the rednced form of the Yang-Baxter 

map (^Urn,m qm,n+l) ^ ?m,n)■ 




-1 


{a* -a*)(l + U„,^nqln,n+l ) - {a - a*) [I + qni,n+iqln,n+l 


1 -1 


X 


® Y ^m,nqm,n+^ ^m,n (® ® Y ?m,n+l?rn,n+l^ ?m,n+l 


qm,n 


(a - a) (J + ^ - (a - a*) (J + Mm.nwL.n) 


-1 


-1 


X 


(cr Cl) (^I + gjn,n+l'^m,n) *?m,n+l (c^ CT ) (^/ + Um,n'^ln,n) 


-1 


U. 


m,n 


Recalling that the parameter a (resp. a) can depend on the discrete inde¬ 
pendent variable n (resp. m), we obtain a nonantonomons system defined on 
the two-dimensional lattice: 

K - Ka) (l + («„+!,„ - 

= (®m ®m) qm,n+iqm,n+lj {'^m+l,n qm,n+l) y 


(cTjj Ot-m') Y qm,n) 

(c^n ®n) ('^ Y {Um,n qm,n) ■ 

This is a fnlly discrete analog of the principal chiral model for Hermitian 
projection matrices in component form: 

= (^ + V) (/ + W*)-' (« - 9). (3.17a) 

= (/ + 5«t) (/ + ( 3 , 17 b) 

That is, fl3.17p is obtained from the principal chiral model fl3.14p by setting 
(cf. (l^ i 

fiv) = -iY(h), ^(0 = 


P 




u 

I 



y 


^ = 





q 

I 
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where a(0, a{r]), ^{r]), e M. 


4 Concluding remarks 

Auto-Backlund transformations for a continuous integrable system provide a 
useful clue to obtaining its proper discretizations, but this still remains to be 
clarihed for two or more component systems such as the NLS system and their 
various reductions. In this paper, following the new approach introduced in 
our previous paper PI. we obtained a new proper space-discretization of the 
matrix NLS system, which admits the Lax-pair representation and permits 
not only the complex conjugation reduction but also the Hermitian conju¬ 
gation reduction between the two dependent variables as in the continuous 
case [13] (see fl2.26p L This is in contrast to the matrix generalization of the 
Ablowitz-Ladik discretization of the NLS system IZ7I, which permits only the 
complex conjugation reduction and not the Hermitian conjugation reduction 
in local form (see [221121] and references therein). Thus, our space-discrete 
matrix NLS system is more appropriate and physically meaningful; indeed, it 
can generate proper space-discrete analogs of various multicomponent NLS 
equations [2111^1 - 136] by considering reductions analogous to the continuous 
case. By changing the time part of the Lax-pair representation, we can also 
obtain proper discretizations of the matrix mKdV system and the matrix 
sine-Gordon system, which admit various interesting reductions such as the 
space-discrete vector mKdV equation 02.2711 . 

In our approach, we reinterpret (a slightly generalized version of) the 
binary Backlund-Darboux transformation as a discrete spectral problem, 
wherein the two intermediate potentials appearing in the Darboux matrix are 
used as a pair of new dependent variables. Bianchi’s permutability theorem 
for Backlund-Darboux transformations implies that the number of discrete 
independent variables can be increased consistently to define a multidimen¬ 
sional lattice. The consistency condition on an elementary quadrilateral can 
be solved explicitly, providing a new Yang-Baxter map of the NLS type, 
which can be regarded as a fully discrete analog of a reduced form of the 
principal chiral model. It would be interesting to investigate whether this 
result can be related to the recent work of Caudrelier and Q. C. Zhang |46] . 
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1997 on multicomponent generalizations of the sine-Gordon equation. 


A Integrable discretization of the vector sine- 
Gordon equation 


In this appendix, we present an integrable discretization of a vector analog 
of the sine-Gordon equation uniEo], which was obtained by the author in 
the late 1990s but left unpublished. 

We start from the continuous case and consider a Lax-pair representation 
of the following form (cf. 02.21) and 02.7P '): 


Ti ■ 


2v^2a; 


■ ■ 

^2 

X 



T 2 


Ti ■ 

i 

1 

1 

2 

Ti 

T 2 


77 

1 

> 

1 

^2 


(A.la) 


(A.lb) 


Here, c is a constant, / is a scalar function and Q and 77 are square matrices. 
The compatibility condition for this overdetermined linear system provides 
the following three equations: 


r/./=(Q77), = (77Q),, 




! -7 2a; 1 4Q — 

VC-/ ) 


Vc- / 


77o 


477 = 


t-i 


O, 

o. 


To satisfy the first equation, we restrict Q and 77 to the form: 

2M-1 2M-1 

Q = g(b/ + 77 = ^ 

i=i j=i 


(A.2) 


(A.3) 


where the x matrices {ei, 62 ,..., e 2 M-i} are generators of the 

Glifford algebra; they are linearly independent and satisfy the anticommuta¬ 
tion relations: 

. GjGk T ‘^^jkl- 


Thus, we have 


Qn = nQ={q,q)I, 
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where q = ..., so the hrst equation in flA.2p is satished by setting 

/ = {^5 q)- Then, flA.2D reduces to the vector sine-Gordon equation 


\/c- {q,q) 


- 4g = 0, 


(A.4) 


t-1 


or equivalently, 


By setting 


^^*-1 + r-In n\ ^y/c - (qr, q) q = 0. 

^ \Q.^ Q) 


-Qx =■ u, 


2^Jc-{q,q) 

which changes the spectral problem flA.laP to a standard form, we obtain 
the potential form of the vector sine-Gordon equation flA.4p : 

1 ,, , 
/■ ; ^=ut_,x = 2w. (A.5) 

^yAc- {ut_,,ut_,) 


More information on multicomponent generalizations of the sine-Gordon 
equation can be found in (TSUBD] and references therein. 

Next, we consider a discrete-time analog of the Lax-pair representation 
flA.lD as given by 


■ 


-Ki 



' ■ 

T2,n . 

X 


KI 


. ^2,n . 


^l,n+l ' 


1 

1 

Qn 

\ 


^2,n-|-l 

c 

T^n 

-y/c- fnl _ 

i 

. ^2,n . 


where Zl is a time-step parameter. The compatibility condition provides 


fn,xl — {QnT^n)x ~ i'^nQn)x ) 

1 

Qn+ly 


^ V Vc - fn +1 


n 


\/c- /r 
1 


■Qn,x 1 — 2Qn+l — 2Qrj — O, 


(A.6) 




-'R-n.x 1 — 2TZn+l ~ 2TZn — O. 


^In constructing explicit smooth solutions, each square-root function in this appendix 
suffers from an intrinsic sign problem. Note that (TOl in the scalar case reduces to 
the sine-Gordon equation by s etting c = 1 and q = sinp, but a sign ambiguity arises in 
extracting the square root as \/ cos^ p = ± cosp. Such a sign ambiguity can be resolved by 
expressing the dependent variables in the Lax pair appropriately in terms of trigonometric 
functions, but then the equations of motion become complicated in the multicomponent 
case. 
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Then, by restricting and to the same form as in flA.3D and setting 
fn = {Qui Qn)j flA.hp simplihes to a proper discrete-time analog of the vector 
sine-Gordon eqnation flA.dp : 


1 

z 


'C- (g„+i,q„+i) 
By setting 

Mr, 


Qn+1,1 


\/c - {Qn. Qn) 


2qn+i - = 0. 

(A.7) 


-Qn,x • '^ni 


2Vc- (qr^,qr„) 

we obtain a time-discretization of the potential vector sine-Gordon eqnation 
flA.Sp as 

'^n,x) '^n+1 '^n) ("^n+l T '^n) ■ 

(A.8) 

Essentially the same eqnation was obtained independently by Balakhnev and 
Meshkov [8Tll8^ as an anto-Backlnnd transformation for a vector analog of 
the mKdV eqnation. 
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